Multiscale modeling and simulation of material phase change problems: ice melting and copper crystallization

Recommended Citation Wei, Xiupeng. "Multiscale modeling and simulation of material phase change problems: ice melting and copper crystallization." MS ABSTRACT The primary objective of this work is to propose a state-of-the-art physics based multiscale modeling framework for simulating material phase change problems. Both ice melting and copper crystallization problems are selected to demonstrate this multiscale modeling and simulation. The computational methods employed in this thesis include: classical molecular dynamics, finite element method, phase-field method, and multiscale (nano/micro coupling) methods. Classical molecular dynamics (MD) is a well-known method to study material behaviors at atomic level. Due to the limit of MD, it is not realistic to provide a complete molecular model for simulations at large length and time scales. Continuum methods, including finite element methods, should be employed in this case. In this thesis, MD is employed to study phase change problems at the nanoscale. In order to study material phase change problems at the microscale, a thermal wave method one-way coupling with the MD and a phase-field method one-way coupling with MD are proposed. The thermal wave method is more accurate than classical thermal diffusion for the study of heat transfer problems especially in crystal based structures. The second model is based on the well-known phase-field method. It is modified to respond to the thermal propagation in the crystal matrix by the thermal wave method, as well as modified to respond to temperature gradients and heat fluxes by employing the Dual-Phase-Lag method. Both methods are coupled with MD to obtain realistic results. It should be noted that MD simulations can be conducted to obtain material/thermal properties for microscopic and/or macroscopic simulations for the purpose of hierarchical/sequential multiscale modeling. These material parameters include thermal conductivity, specific heat, latent heat, and relaxation time. Other type of interfacial parameters that occur during the phase change process, such as nucleus shape, iv interfacial energy, interfacial thickness, etc., are also obtained by MD simulation since these have so far been too difficult to measure experimentally. I consider two common phase change phenomena, ice melting and copper crystallization, in this thesis. For the case of ice melting, MD is first employed to study its phase change process and obtain thermal properties of ice and water. Several potential models are used. I conduct simulations of both bulk ice and ice/water contacting cases. It is found that various potential models result in similar melting phenomena, especially melting speed. …